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We introduce the quantum Levy walk to study transport and decoherence in a quantum random 
model. We have derived from second order perturbation theory the quantum master equation for 
a Levy-like particle that moves along a lattice through hopping scale-free while interacting with a 
thermal bath of oscillators. The general evolution of the quantum Levy particle has been solved for 
different preparations of the system. We examine the evolution of the quantum purity, the localized 
correlation, and the probability to be in a lattice site, all them leading to important conclusions 
concerning quantum irreversibility and decoherence features. We prove that the quantum thermal 
mean-square displacement is finite under a constraint that is different when compared to the classical 
Weierstrass random walk. We prove that when the mean-square displacement is infinite the density 
of state has a complex null-set inside the Brillouin zone. We show the existence of a critical behavior 
in the continuous eigenenergy which is related to its non-differentiability and self-affine characteris- 
tics. In general our approach allows to study analytically quantum fiuctuations and decoherence in 
a long-range hopping model. 
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Random walks have been studied for a long time now in order to obtain information about the effect of dimensional- 
ity, symmetry, and topological structures on the general properties of classical transport phenomena. Via the Central 
Limit Theorem the Gaussian distribution plays a fundamental role for all random walks (dynamic semigroup) with 
QQ finite mean square displacement per step [l|, y| . Levy fiights are classical Markovian random walks in a continuous 
space with infinite mean square displacement per step (thus Levy fiights are not in the attractor of the Gaussian 
distribution). The interesting point about these walks is that the set of points visited by these Levy flights has a 
self-similar clustering property, so a fractal dimension has been extensively discussed and in particular magnificently 
' illustrated by Mandelbrot 3]. A Weierstrass random walk is a discrete (lattice) version of Levy fiights which was 
T-H ■ introduced by Schlesinger et al. [1]. This transport model (also called Levy walk) shows the fractal nature of the 
^ ■ random walk trajectory, and several mathematical concepts such as non-differentiable function, lacunary Taylor series, 
. i-H , fractals, renormalization group transformation, etc., have been related to this Weierstrass random walk and mapped 
to important characteristics of transport phenomena fsj . 
Jh Recently it has been reported that power-law discrete-stable (step) distributions can be obtained as the stationary 

. . state of immigration Markov processes [6| , showing that the statistics of the steps of the particle is inherited from the 
TO— tuplets rates (mesoscopic) fiuctuations appearing in the master equation [7]. To tackle this kind of program from 
a non-classical point of view is much more harder because in quantum mechanics it is necessary to include the cause of 
damping and noise explicitly, thus apart from the Hamiltonian of the system of interest S one has to include a thermal 
bath B an interaction between both P, Q . The program of the present work is not in the direction of reference , 
rather our purpose here can be presented in two fold: first to get a quantum semigroup description of a mechanical 
particle that moves doing scale-free hopping along a lattice of sites, and second to study the quantum decoherence 
phenomena in a long-range dissipative model. In order to achieve this program we will adopt the simplest quantum 
mechanical assumptions to get an analytical model. 

To gain insight into the relationship between dimensionality, topological structures and disorder, many quantum 
models of transport phenomena have been introduced [9j. Among them the tight-binding approximation for a quantum 
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particle over a regular structure with nearest-neighbor (NN) interactions, is a simple description which is equivalent 
to the quantum random walk (QRW), a quantum particle that moves along a lattice of sites doing NN steps while 
interacting with a bath [1, ■ One of the most interesting facts that distinguishes quantum mechanics from classical 
mechanics is the coherent superposition of distinct physical states. Many of the non-intuitive aspects of the quantum 
theory of matter can be traced to the coherent superposition feature. Two important question are: 

• How does the coherent superposition operate in the presence of dissipation? 

• How does a long-range interaction characterize quantum decoherence? 

These subjects have been important issues of research since the pioneer works of Feynmann and Vernon (ill, [l^. 
Caldeira and Leggett [l^ among others, see for example the references cited in: [ll. [l3 - [l7| . 

The study of a quantum walk subjected to different sources of decoherence is an active topic that has been considered 
by several authors, in particular, by their interest in understanding Laser cooling experiments p^ . modeling Blinking 
Statistics [l^, and also in doing quantum simulations [l^l- It should be noted that the most usual Laser cooling 
scheme is based on the idea that the microscopic quantum description of subrecoil cooling can be replaced by a study 
of a related random (anomalous) walk in momentum space, this is the point where the concept of Levy waiting-time, 
for the statistical description of the elapsed time between walks, appears. The definition of our Quantum Levy Walk 
(QLW) is base in a long-range jump model without introducing any waiting-time statistics. The physical motivation 
of using a Levy-like probability for the jump was inspired on recent experiments showing long-range interaction as in 
Rydberg gases |2l| - [23| . We end this paragraph noting that our QLW is quite different from the quantum walk with 
Levy waiting-time. 

In this paper we introduce a scale free open quantum model, i.e., a quantum mechanical particle that moves along 
a lattice through hopping scale free while interacting with a thermal phonon bath. We have chosen the interaction 
Hamiltonian with the bath in such a way that it produces a long-range superposition of vector states. We highlight 
some of the issues of interpretation of the coherent superposition by tackling a soluble long-range hopping model. 
The asymptotic long-time regime of the quantum purity is characterized by a long-tail with a non-trivial exponent 
that depends on the Weierstrass parameters of the Hamiltonian. A long-time coherent behavior for the localized 
correlation function is also explained in terms of the present scale free hopping model. 

In appendix A we present some formal aspects of quantum dynamic semigroups, and we revisit the second order 
approximation obtained for an open quantum system weakly coupled to the environment [l|, Il^ - [l7l [23 - l27| . Thus we 
emphasize some general conditions on the system of interest, the environment and the interaction Hamiltonian to 
obtain a true semigroup. In the core of the paper we present our open quantum model: The Quantum Levy Walk, 
which is a generalization of the QRW Hamiltoniaii [Ml ■ Then we obtain the evolution equation for the reduced density 
matrix under the Markovian approximation (28l . [29| . its solution, and solve analytically some correlation functions 
associated to the coherent superposition feature. 



For open quantum systems, the Markovian description of the dynamics is based on the concept of quantum dynamic 
semigroups [l^ [IJ] . Only with these semigroups are the properties of the reduced density matrix of the quantum 
system of interest, preserved during the whole time evolution (positivity, trace and hermiticity) (25l - [27j . From a 
microscopic description considering the total Hamiltonian of the system of interest and the environment, it is possible 
to derive a picture involving the quantum dynamic semigroups. 

One of the pioneers work in obtaining the QRW model from first principles can be found in van Kampen's paper , 
where it is shown that tracing out the baths variables a bona-fide semigroup is obtained. Here we will do something 
similar, but introducing a scale free Hamiltonian, for the free particle, and generalizing the interaction with the phonon 
thermal bath. Let the system 5 be a free particle that can reside in any lattice site / = 0, ±1, ±2 • • • , the dynamics 
of the system S will be described by the Hamiltonian 



where the shift operators a, act on the orthonormal set of Wannier basis |/) that spans the Hilbert space H. A 
general form for these shift operators can be expressed as: 



II. THE QUANTUM LEVY MODEL 
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with /(e„) > and fulfilling normalization to one, to help its physical interpretation. 

A particular choice of /(cn) will describe a tight-binding-like Hamiltonian ranging form NN to long-range hopping, 
thus the set a, will be the fundamental operators to modelling the interaction Hamiltonian between system S and 
the bath (a thermal set of oscillators), see appendix A. A classical one dimensional random walk is defined in terms 
of the probability for a particle to make a step of a given length to the left or to the right. Quantum random walk 
Hamiltonians are described instead in terms of probabilities amplitudes (here the shift operators a and a^). Related 
NN discrete-time models are (coined) quantum random walks [30| . Also the NN random walk Hamiltonian (in a 
ring) has been used to study transport in a quantum trapping model (3l| . Quantum random walk, is the counterpart 
of a classical random walk for particles which cannot be precisely localized due to quantum uncertainties. When 
/(fn) — 0, Ve„ 7^ €o with lattice parameter eq = 1 we reobtain the usual NN random walk in the line Q. When 
/(e„) oc n^^^", en — n, a > 0, n — 1,2,3 ■ ■ ■ we get the Gillis and Weiss lattice walk model [111, etc. 

If f{tn) is characterized by a power-law probability we obtain a Levy-like jumping walk, this type of quantum 
walk has been previously reported in order to study long-range interaction in a non-dissipative model (ssj . as well 
as in a trapping transport model [s^l- Here we propose to study a jumping model characterized by the Weierstrass 
probability, this class of lacunary long-range probability has been widely study in classical transport phenomena, 
nevertheless up to our knowledge nothing has been reported in the context of quantum transport. Two important 
consequence appear from the definitions ([2]) and ([3]): using that e„ belongs to the Wannier index we can see that a 
and are diagonal in the Fourier basis, this fact leads to the conclusion that also Hs will be diagonal, and also it 
can be proved that a and commute. These results will be shown in the Appendix B for the particular Weierstrass 
jumping probability. 

We define the Weierstrass shift operators in the form: 

A-l °° °° 



E E i 1^ - ; > 1' ^integer) > 1, (4) 
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«^ = E E 7^ 1^ + (^1 ; > 1' ^integer) > 1. 
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This means that the application in the Weierstrass shift operators to any vector produces a linear combination 
in the Wannier basis. This linear combination is scale free and has a clustering structure [STj . This clustering is 
characterized -in average- by the probability cx 1/^" to have a projection on a Wannier vector distant 6" from 
For example 
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Interestingly as we mention before. 



[a,at] = 0; Vyl > 1, V6(integer) > 1, 

this and other important results are shown in Appendix B. 

We noted that only in the case 6=1 the operators a and arc truly translation operators, in the sense that 
the product of successive translations is equivalent to one resultant translation, i.e., for example the application of n 
translations gives: 

a'^ ■■■a)\l) ^\l^n) only if 6=1. 

Therefore it is simple to see that taking 6 = 1 in ([4]) we reobtain the QRW model [1, [13] • 

In general for 6 > 1 the Hamiltonian Hs describes a scale free tight-binding-like Hamiltonian (the QLW model), 
i.e., a quantum free particle in a lattice that moves, making hopping like a Levy walk, while interacting with a bath 
(in the present paper we work in a discrete one dimensional infinite Hilbert space). The eigenfunctions of Hs are 
denoted by the kets \k) (with — tt < A; < tt) and are given by the Fourier transform of the Wannier vectors: 



oo 



OO 



4 



thus (fell Hs Ife) = £ki {b, A) S (fci — fc2), where 

c ^^ r^ /i -^-^ COS (6"fc) 1 , , , 

gfc(6,^) = 0|l--^2^^ ^ 1 ;-^<fc<7r. (5) 

The eigenenergy ^) can be related to lacunary Taylor and Fourier series |35i|. From ([5]) we can define A(fc) = 
1 — j^£k{b, A), this function obeys the scaling equation 

A(fc) = ^X{bk) + cos(fc). 

It has been shown 36] that the nonanalytic part As (A;) of A(A:) satisfies the homogeneous equation 

As(fc) = ^A,(6fc), (6) 

so that 

\{k) = \k\^ Q{k), with /i = Inyl/lnfe, 

here (5(fc) = Q{bk) is a bounded function periodic in In |fc| with period In 6 of quite intrincated structure Therefore 
our eigenenergy £k {b, A) shares some similarities with critical phenomena analysis (renormalization group transfor- 
mation). In addition if the parameter b is an integer it has been shown [35] that the power series in 2: = e**^ of A(fc), 
has gaps or missing terms. These gaps lead, by using the Fabry's theorem, to the concept of noncontinuability of the 
series of X{k), and so to the conclusion that A(fc) has extremely complicated behavior as a function of k. In fact, for 
/i < 1, A(fc) is Weierstrass' example of a function which at no point possesses a finite derivative. When this result is 
translated to the eigenenergy £k{b,A) we may conclude that the density of states (DOS) is not well defined if /i < 1 
this will be shown also in the next subsection, see Eq. ([5]). 

As we mentioned before the type of Hamiltonians ([Ij, with shift operators as presented in ^ and share the 
properties of been diagonalized in Fourier space. The important point of our Weierstrass quantum model is that the 
eigenenergy turns to be non differentiable for the critical value fj, < 1 (i.e., b > A), this fact will be analyzed in the 
context of the time evolution of the reduced density matrix, in detail, in the next sections. 

An analysis concerning quantum walks with long-range steps (but without dissipation) was carried out by Miilken 
et al. (33| , in that paper it was shown that there exist a universal behavior for the quantum walks which is different 
from the universality of long-range classical random walks. In our present work we go one step forward and study 
a long-range quantum walk coupled to a thermal bath, then new universalities are found for the decoherence of the 
system. 



1. Density of states 

Using the Green function of the Hamiltonian Hs the DOS can straightforwardly be calculated if we known the 
matrix elements of G{Z) = [Z — Hs]~^ , for example by using the formula D{£) = \ims^o+ -^^TxiTr [G{£ + iS)] . 
Nevertheless for the particular case of our QLW Hamiltonian, and due to the fact that we already have an analytic 
expression for the continuous energy eigenvalue £k{b,A), it is more convenient here to calculate the DOS by the 
alternative formula 



D{£) 



1 

2^ 



S{£~£k(b,A))dk 



(7) 
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where kj = kj{£) are the solutions of the transcendental equation: £ = £kj{b,A). From ([5]) we see that 



d£kib,A) ^^A 



dk 



A 



-J 00 



n=0 



sin (&"fc) 
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then we expect that when b > A the right-hand-side may diverge for some values of kj. This fact ultimately leads 
the DOS to be not well defined, this issue can be seen as the occurrence of a complex null-set inside the Brillouin 
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FIG. 1: Plot of £fc(fe,^)/r2, the continuous energy eigenvalue of Hamiltonian H]) as a function of the Fourier number k in the 
first Brillouin zone (reciprocal to the Wannier lattice index). In (a) and (b) the straight lines correspond to the tight-binding 
case, i.e., the QRW model (6 = 1). The other lines (dotted, slashed and dotted-slashed) correspond to the QLW for different 
values of Weierstrass' parameters: b{= 2, 4), ^(= 2, 4, 6). For large values of the Weierstrass' rate b/A > 1 the non-differentiable 
structure of £k{b,A) is clearly visible as was predicted from the scaling ([6]) and also from Eq. ((8]). 



zone, and is in agreement with the previous mathematical report, coming from ([5]), on the non differentiability of A(fc) 
when n < 1. Only when b < A the eigenenergy £kib,A) is differentiable anywhere and the behavior of the DOS for 
the Weierstrass' model shares analogies with the DOS numerically calculated by Miilken et al. in a long-range model 

The Hamiltonian ([T]) in the limit 6 — > 1+ is the QRW model [1, [13, HH, then as expected, we can reobtain from 
([5]) the usual tight-binding eigenvalues £k{b = 1,A) = ~ cos (fc)}, and from ([7]) the corresponding DOS: 

D{£) ^ ^ {2ns - £^y^^^ , 0<£<2n if b = l. 

In conclusion: for 6 > 1 the DOS follows from ([7]), but this density only is well defined under the constraint b/A < 1. 
In the opposite case b/A> 1 the DOS is not defined because the function £k{b, A) is non-differentiable. In figure 1 we 
have plotted £k{b,A) for several values of Weierstrass' parameters b,A, in this plot the non differentiable signature 
of the eigenenergy of the Levy walk Hamiltonian can clearly be seen. 

Here it is important to note the difference with the classical characteristic function of the Weierstrass random walk 
(discrete space Levy flight 04Sls3])- the classical walk the condition to have an infinite mean square displacement 
per step is b^/A > 1 (because the second moment is given by the second derivative of A(fc)). Then when clustering 
occurs, the number of subclusters within a distance &" of the origin is, on average, A'"'. Note that for the classical 
walk, the characteristic function of the Weierstrass walk is non-differentiable when < /i < 1. Then it is possible 
to call to the value /i = In A/ In b the Hausdorff-Besicovitch dimension of the set of sites visited by a classical walk 

35I [37I . In quantum mechanics a similar analysis could be done in the context of the Schodringer-Langevin picture 

28| . i.e., without using the density matrix. 
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We want to remark that if /i < 1 (6 > A) the eigenenergy £k{b,A) is non-differentiable and the DOS is not well 
defined. A rigorous calculus based on the scaling properties ^ of the eigenenergy leads to the conclusion that the 
record of £k{b,A) as a function of k, for b above the critical value b > A, is a. self-affine function. In fact, a fractal 
dimension can be measured (for example) by using the Box-counting technique and the prediction gives D ~ 2~ fi for 
< ^ < 1 [3, m, [3^ . An alternative technique based in the analysis properties of the zero crossing 4Ql has recently 
been reported to be suitable to characterizes processes that are continuous, but the derivative has fractal properties. 
Interestingly many years ago these processes were termed "subfractals" I41'|. 

Note that we are assuming integer values for b to avoid non-commensurability problems for the walk in the associated 
discrete infinite dimensional Hilbert space H of lattice parameter e = 1 (at the end of appendix B we present some 
discussion on the case when the lattice parameter e goes to the continuous limit). 



A. The Quantum Master Equation for the quantum Levy walk 

Concerning the quantum bath (see appendix A for details), we will assume that the thermal bath B is an infinite 
set of oscillators, and that the interaction with the phonon bath causes the free particle to jump either to the right or 
to the left. Thus we consider the operators Vg appearing in the interaction Hamiltonian (IA8[) to be proportional to 
the Weierstrass shift operators a and defined in Eq. In the present paper we will study a particular coupling 
with the heat bath, nevertheless going back to appendices A, B and C it is straightforward to write down the QME 
considering any other type of coupling. For example: if Vg oc a — we can model an interaction that in the continues 
limit (lattice parameter going to zero) coupled the velocity of the free particle with the heat bath. 

Taking into account that [a\ a] — we trivially get that in the Heisenberg representation Weierstrass's shift operator 
does not have a time evolution, i.e.,: 

Vpi-T) = eM~^rHs/h)V^eM^rHs/h) = VfsiO). 

This result ultimately will lead to the fact that the Kossakowski-Lindblad (KL) generator [25|, [2^ will be completely 
positive [2^. Thus, using Eqs. (jA9p and (|A10I) in (|A3I) . we can write the Quantum Master Equation (QME), for the 
reduced density matrix p, in the form (see (|C9|) in appendix C) 



p = — ^ [Heff,p] + ^ {2apa^ - a^ap - pa'' a) + ^ (2a Va - aa^ P - paa^) 



A(3h 



(9) 



Here a > is the dissipative constant and /3 the inverse of the temperature of the bath. 

Note that in the case 6 7^ 1 we have a' a 7^ 1 (see appendix B) therefore the effective Hamiltonian has a non-trivial 
contribution: 



H^ff — Hs — hjjco'a^ (jJc > 0. 
The effective Hamiltonian can also be diagonalized in the Fourier basis, therefore we get 



(10) 



Ik.) = 



£k^ (6, A) — huj. 
EkAb,A) S{ki~k2 



A-1 
A 



E 

n2,ni—0 



cos fci (6"=^ - 6"! ) 



S{ki 



(11) 



where we have used (jB2p and E^-^ {b, A) is given in ([5]) . Note that HlUc is an upper bound energy which is characteristic 
of the Ohmic approximation [13], (see (|C8p in the appendix C). In figure 2 we have plotted Efe(&, A) for several values 
of Weierstrass' parameters b, A, in this plot it can be seen that the effective eigenenergy associated to the Hamiltonian 
PH)) has larger fiuctuations than E^^ {b, A) the eigenenergy associated to the naked Levy walk Hamiltonian. 

When b — 1 (the QRW case) the effective Hamiltonian only introduces a trivial constant because limb_j,i a^a 1. 
Therefore, for b — 1 the KL generator has a very simple expression, form © the QME results [1, [l^] 



L[p] 



[Hs,p] 



2(3h 



{apa) 



a' pa) 



if 



the group iTa/2(3h can be called the diffusion constant D and is given in units of [time] ^ because the lattice parameter 
is dimensionless. The case _D — 00 corresponds to the high temperature limit, see appendix C. 
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FIG. 2: Plot of Eki{b, A) /fl, the continuous energy eigenvalue of the effective Hamiltonian (|10p as a function of the Fourier 
number k in the first Brillouin zone. In (a) the energy rate is re = hujc/^l — 1, and in (b) is re = 2. The straight lines 
correspond to the tight-binding case (6 = 1). Other lines (dotted, slashed and dotted-slashed) correspond to the QLW model 
for different values of Weierstrass' parameters b{= 2,i),A{= 2,4). The physical meaning of the contribution hujca^a in the 
effective Hamiltonian H^ff can be thought as increasing the multiplicity of the hopping structure, see (|B1|I . For b > A and 
large values of the rate the amplitude of the non-differential structure get larger. 



B. On the second moment of the quantum Levy walk 



From the QME ^ we can obtain the dynamics of any operator, in particular here we are interested in the evolution 
of the dispersion of the position operator q, which in the Wannier basis has the matrix elements: 

{h\q\h) ^lSi„i,, (12) 

note that q is defined as a dimensionless position operator. In the \l) representation it is possible to see that the 
thermal mean-value time-evolution of the first and second quantum moments can be written in the form: 

|(q(t)) = |rr{q(t)p(0)} = rr{qp(t)} (13) 

= ^rr{p[q,iJe//]} 

^(q\t)) = rr{qV(t)} (14) 



dt 



^ ' Til ,n2=0 ; = — CXD 
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Setting 6 = 1 in (|13p and (|T4)) . i.e., in the limit of the usual QRW, we can write in the Heisenberg representation: 

These equations can easily be solved. Using that [q, Hs] — ^ {o- ^ a^) (see appendix B), and that a{t) — a{0),a^{t) ~ 
a^(0) we get for the time evolution of the position operator: 

q{t) = ^{a-a^)t + q{0) if 6=1. 

In the same way, the variance of the QRW (the NN walk model) can be calculated giving: (q(i)^) — (q(t))'^ — 

1 + 2Dt, with 2D = ^, which is the expected dissipative resuh [H, i, B [H . From ^ and (US]) it is possible 

to see that von Neumann's term gives a contribution of the form cx for the time-evolution of the second moment, 
this is a well known quantum result [sl], see eq. (pUj) for a general discussion. The dissipative contribution comes 
from the interaction with the bath B, giving the classical diffusive behavior cx t. 

From (1141) and due to the coherent dynamics involved through the time-evolution of the off-diagonal elements 
{l \ p{t) \l — -I- 6"^), it is not simple to realize what will be the dynamics in the general case when b ^ 1. But in 
principle any higher moments of q(t) can also be analyzed in the same way from our QME ([9]). 

Let us now analyze the case 6 7^ 1. In this case the interaction with B produces long-range hopping, and consistently 
a non-trivial quantum decoherence phenomenon. The dissipative term of (jl4p . can also be written in the form: 

i ^ E ^ E « + .»-i.wi<.."->. (17, 

^ ' Tii,n2=0 ^ ' /= — 00 

We can explicitly calculate this contribution by going to the Fourier representation. First of all, here we will assume 
that the initial condition for the reduced density matrix was prepared in a pure Wannier state: p(0) = |^o) ('o|i so 

piO)k,M = (fci|p(0)|fc2> = ^exp{-i{ki-k2)lo): (18) 

thus {k\ p(0) |fc) = On the other hand, from ([9]) is possible to see that ^ (fc| p(t) \k) = therefore (fc| p(t) \k) = 
^,Vt, (see next section). Now going back to ([T7|) we can write: 

/ A—^ \ ^ r r °° / h \ "^^"^ qq 

m-^i^) dk.dk, [a] E (^+^"i^i)(fcii/^wi^2)(fc2 1^+6-) 

Trn / A-}\^ r r ^ //)\"'^"' ^ pi{ki-k2)l 

TV(X / J\ — i \ II,,,, X ^ / f \ _if,"2 1,„ ,, I , I, , X ^ c 



pnyA J J J ''^''^ ^ 1:4) Mfc,|p(.)|fc.) E 2. 



ni ,712—0 



ni,n2— ^ 



noting that we have assumed h = integer we can write: 



_ TTo M-iy / ^ sin(b"i -6"^)7r 

— if 6=1 

ft^^, ^> l,6(integer) < A (19) 

00, if b > A. ■ 
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The expression (|19p shows that it is only in the case when b > A that a divergent behavior for the thermal second 
moment of the QLW, (q^(i)) , can arise. This quantum result is quite different from the classical clustering Levy 
flight counterpart (6^ > A) 0-01131. 

The explicit solution of (q^(i)) can alternatively be obtained calculating Tr [q^p(t)] (this is done in appendix D): 



1 fn. 



A^t) +JEt + ll -4 > 1,6=1, 

(q'W) = <! (20) 

''\kw*' + TTri7^l0rW*' + jK')+il A>l,b<A. 



(A-iy 



By taking the limit of null dissipation, i.e., a = we reobtain the quantum behavior {^q^{t)') cx for b < A, this 
result is in agreement with previous reports on the universal behavior of quantum walks with long-range steps |33| . 

We conclude this section noting that the quantum coherence enlarges the threshold (compared to the classical one) 
to have a finite second moment for the walk. For example, if & = 2,^ = 3 the classical second moment is not defined, 
but the thermal quantum average (q2(i)) is finite! 



C. Time evolution of the density matrix 



1. Off-diagonal relaxation 



The QME ([9]) can be solved in the basis of the eigenvectors of Hs- Introducing the representation \k), and using 
that (fell Heff \k2) = Efcj S (fci — k2), with E^^ = F,ki{b,A) given in PT|) we get (see appendix C): 



{ki\p\k2) = — (Efe, -Efej(fci|p|/C2) 



TTa ( A—\ 



2l3h V A 



2 oo 



E :4;r^[2eos(fci6--M-) 

7li ,712—0 

- cos/ci (&"i - 6"^) - cosfca (&"' - 6"^)] (fci| pl^s) • 



From ((2T|) the general solution for the reduced density matrix is 

{ki\p{t) \k2) = p{0)kik2 exp{J^ {ki,k2,b,A)t) , 



with 



T (fci, ^2, b,A) 



E ITT (^"^i) - (^"^2)) 



A 



A-1 
A 



na ( A—\ 



V A 



E 

n2,ni— 
2 oo 



COS kx (&"=^ - 6"! ) COS fca (6"" - 6"^ ) 



E 



1 



[2cos(fci6"i - fc2 6"=) 



- cosfci (6"^ - 6"^) - cos/ca (6"^ - 6"')] , 

note that T (fc, fc, 6, ^) = 0. 

For the QRW (case 6=1) the QME ((2T|) reduces to a simpler form 



(fell P 1^2) = 



- — 1 7T(y, 

— (Efc, - EfcJ + — (cos(fci - k2) - 1) 



(fell P 1^2 



(21) 
(22) 



(23) 



where E^ = ri{l — cosfc} — hujc- Note that in the NN case, the off-diagonal relaxation of the element Pk,o{t) (from 
the homogeneous Fourier mode) is just controlled by a trivial trigonometric function: 



T{k,0,b^l,A) = - 



*7i ^ ^ 



(1 — cos/c) 
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This result resembles the relaxation of the Fourier modes of a classical one dimensional ordered walk with a diffusion 
coefficient 2D — ^ (see appendix E). Note however that in quantum mechanic the relaxation of the initial position 
is controlled by a double Fourier integral, so we cannot expect the same long-time asymptotic behavior as in classic 
(see next sections for details). 



2. Diagonal relaxation 



In general for & > 1, by putting fci = ^2 in (|21l) . we see that the diagonal elements pkk{t) remain constant in time: 

l{k\p\k)^0. (24) 

This result tell that even when there is a decoherence in the off-diagonal elements of the density matrix, the probability 
distribution of the Fourier modes is invariant in time. The physical interpretation of this fact can be understood by 
noting that the thermal average of the kinetic energy is constant in time. As we mention before this is a consequence 
of the general form of the shift operators ([5]) and 

To prove this fact for the Weierstrass model, we first define a pseudo-momentum operator: 

TT) 

P^-7:[ci,Hs], (25) 
in 

where m represents the mass of our free particle in the lattice. From (jB3[) it is simple to see that p is diagonal in the 
Fourier basis p |fc) = pfc |fc), i.e., 

{ki\p\k) = !!^(^:l^^^(^Aysin(fc6")<5(A:-fci) (26) 
= S{k~ki)pk. 

Once again the eigenvalue p^ is well defined only below the threshold b < A. 

Assuming the system was prepared in the pure state: p(0) = \Iq) {Iq\, the thermal average of the pseudo-momentum 
gives (p(i)) — 0,Vt, and for (p^(<)) we get 

(p^t)) ^ Tr[p^p{t)]^ f dk f dki{k\p^\ki){ki\p{t)\k) (27) 

J —TT J —TT 

dk f dki {pkfS(k,-k){ki\p{t)\k) 

1 

= — / dk (ph) = constant, 
2tt ^^^^ 

where we have used that {k\ p{t) \ k) = -^yt. In general it is possible to prove that the thermal average of any 
observable which is diagonal in the Fourier basis, will be constant in time. 

It is interesting to comment here that for the QRW (the NN case) the pseudo-momentum operator (|25p is in fact 
a discrete version of the momentum, and can be written in the form (see appendix B) 

P"7ft2(""")' ^=^' 

then it is simple to see that p = ^ [p, Hs] — 0, so the discrete momentum in the Heisenberg representation is constant 
in time. The "lattice" commutation relation between the position and the momentum gives 

[q,p]^-—{a + a^), if 5=1. 

Thus in the continuous limit (taking the lattice parameter e — >■ 0, see appendix B) we reobtain the usual commutation 
relation [q, p] = ihl, p^ . 
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3. Quantum decoherence from a pure state 

In order to study the quantum decoherence due to the interaction with the thermal bath B, we propose here to 
analyze the time evolution of the density matrix assuming that at time t = the system was prepared in a pure state: 
p(0) = |/o) (^ol- Then the quantum probability to be at site I at time t is given by 



{l\p{t)\l) = 



dki 



2 



dk2 {I\ki){ki\p{t)\k2){k2\l) 

dk2 exp [i {ki - ^2) {h \ pit) |/c2) 
dki / dk2exp[i{ki — k2) {I ~ lo)]exp{J^ {ki,k2,b,A)t) . 



(28) 



1 

27r, „ „ 

Calling (/ — Iq) = Al the distance form the initial condition, we can plot the probability P{Al,t) = {l\p{t) \l) as a 
function of Al for different values of time t, and Weierstrass' parameters b,A. Note that P{Al,t) < 1 as expected 
because the density matrix is well normalized (^1 p{t) \l) = l,yt ^ 00, as can easily be checked from ((28)) . In figure 
3 we plot P{Al,t) for different values of the frequency rates r = j;/D, = uJc/D and Weierstrass' parameters b,A^ 
the non diffusive characteristics of the profile P{Al,t) as well as its trimodality behavior can clearly be seen. This 
trimodality is the result of the combination between the non-diffusive characteristics (clustering) of the Levy walk 
(by increasing b > A) and the quantum oscillations (coherence by increasing r > l,rc > 1). All Fourier integrals were 
made using the numerical integration quadrature method 42] . 

From (j28p we see that in quantum mechanic the relaxation of the initial position is controlled by a double integral: 



{lo\p(t)\lo) 



dki / dk2 exp [J- {ki,k2,b, A) t) 

-TT J —TT 



For the case b — 1, i.e., the QRW, we can calculate analytically the long-time asymptotic behavior in the following 

way 



{lo\p{t)\lo) 



dki 



dki 



dk2 exp 
dk2 exp 



h 
in 

li 



(cos ki 



EfcJ + ^(cos(fci 



(29) 



cos fc2) 



(cos (fci 



^ _l51t cos fc^gj^p 



^ (C0s(fci - k2) - 1) 



t dfco. 



it is 



In the limit of t — > cx) we can use the method of steepest descent [43|, so defining the parameter r = ^/ 

possible to see that for r 7^ each Fourier integral introduces a factor l/Vrt, thus the relaxation (in one dimension) 
goes like oc \/rt. On the other hand, if r = (infinite temperature limit) it is also possible to see analytically that 
the dominant contribution in the double integral comes from a small area near ki = fc2 , then we get for the quantum 
asymptotic behavior: oc l/^/i, like in a classic random walk (see appendix E). 

In general for 6^1 another alternative to study the quantum decoherence, is to evaluate some correlation function 
associated to the interference measurement phenomena. This can be done by analyzing the localization probability 
\ p{t) From the general expression for the density matrix, p{t), in an arbitrary basis \4>j), we have 



p(t) = ^a(i)q(t) 10.) (^,1 



then we may conclude that \Cj{t)\^ \{4'j \ ^o)|^ represents the classical probability mixture to measure the localiza- 
tion at time t. We can take into account the relaxation of the quantum interference, from the pure state p(0) — \lo) 
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FIG. 3: Plot of the position probability P{Al,t) when the system is prepared in the pure state p(0) = \lo) {lo\ as a function of 
Al (dimensionless space difference Al = I — lo), for several values of dimensionless time t {D is given in units of [time]"^). In 
the insets we show the QRW case (6 = 1), then it is possible to see that by increasing the frequency rate r = ^/D the quantum 
time-oscillations around AZ = persist longer before the profile spreads in a diffusive way. In figures (a)-(d) we show the same 
profile but for the QLW model (fe 7^ 1) for r — ^,1, two values of the Weierstrass' rate h/ A and frequency rate Tc — u)c/ D. 
For the QLW there is a non-trivial effective Hamiltonian (jlOl) . then we have to specify Tc which is the dimensionless number 
comparing the Caldeira-Legget upper frequency against the diffusion coefficient. In (a), (d) and (c), we chose b/A—l then 
the trimodality of the profile is easily seen. In (a) for Tc = \ the non-diffusive behavior is more visible than in (c) for = 1, 
this trimodality behavior is enhanced by increasing the Weierstrass' rate h/ A > 1 and the frequency rate r. In (b) we chose 
b/A < 1 so the second moment of the P{Al,t) is finite, see (|20p . here the profile spreads but its non-diffusion characteristics is 
clearly visible. 



(ZqI, as follows: the probability to measure at time t the particle at the initial position is from 



{h\p{t)\lo) = \^] I dh I dk2[l-S{h-k2) + S{h-k2)]c^p{T{h,k2,b,A)t) 

1 

2^ 



dki / dk2[l — S {ki — k2)]exp{T {ki,k2,b,A)t) 

J — 7r 

dki exp {T (fci, fci, 6, A) t) . 

Thus using that exp (J^ (fci, fci, 6, ^) t) = 1, we see that the quantum interference is characterized by the localized 
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FIG. 4: Plot of the localized correlation xi^) for the QLW as a function of dimensionless time t when the system is prepared 
in the pure state p(0) — \lo) (lo\. The function x(i) is plotted with r = i for two values of Weierstrass' rate b/A and frequency 
rate rc. The insets show the log-log asymptotic behavior of Xo(i) = x(*) + ^ against the QRW decay ocl/t (wide straight 
line), also in the inset we have plotted (thin straight line) the classical RW localized probability Po(i) oc 1/y/t, see HElf) . The 
long-time coherent behavior of x(i) foi' the QLW is much more visible for large values of rc. Note that if r 7^ none of the 
cases: QLW or the QRW behave asymptotically like a classical ID walk (see appendix E for details). The localized quantum 
correlation x{t) is bounded from above oc 1/y/t and its wavy behavior for long-time (shown in the inset) is the signature of the 
slow decoherence due to Levy's hopping structure. 



correlation function x{t) — (^o| p{t) l^o) ~ i-e-, 

X(i)=(^^^ / '^^^ j ^fc2exp(J-(fci,fc2,&,^)t)-^. (30) 

In figure 4 we have plotted the localized correlation x{t) a function of time for different values of Weierstrass' 
parameters 5, A, and for several frequency rates r, Tc characterizing different energy regimes in the QME. In this 
figure the long-time coherence of %(<:) for the QLW can be compared against the diffusive decoherence corresponding 
to the QRW case, i.e., o^l/t. On the other hand, only in the infinite temperature limit r — the classical behavior 
oc l/^/t is reobtained (see Appendix E). There is some numerical evidence that x{i) for the QLW is in fact bounded 
from above cx 1/ -v/i , and its long-time wavy behavior is the signature of the slow decoherence due to the scale- free 
hopping structure. 

In order to study the decoherence of a pure state due to the interaction with the thermal bath _B, we propose here 
to analyze the time evolution of Tr [p^] , this function is some time called the quantum purity. Assuming that at time 
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t = the system was prepared in the Wannier pure state p(0) = \Iq) {lo\, the quantum purity is given by 

Tr[p{t)^] = f dh f dk2{ki\p{t)\k2){k2\p{t)\h) 

J —7T J —TZ 

= (^27r) J '^^^ J '^'^2exp[J"(fci,fc2,6, yl) + expJ^{k2,ki,b,A)]t. 
Introducing the definition of the rates r, rc in the cxpUcit expression of F (fci, k2, b, A), from (|23p we can see that 

Tr[p{tf] = (J-^ J dki j dk2 ex^v{F{k^MAA)l^^^^^2Dt]. (31) 

Thus the time can be rescaled in unit of the diffusion constant -D, we also see that the behavior of the purity will be 
similar to the localized correlation function x{i) but subtracting the quantum time-dependent oscillations, this is so 
because F (fci, ^2, 6, .4) is evaluated at r = rc = 0, compare with (PO]) . In figure 5 we show Tr in a log-log plot 

for different values of Weierstrass' parameters 6, A. From this plot the non-trivial long tail behavior of the QLW can 
clearly be compared against the QRW case. For the case 6=1 (the NN walks) and using the method of the steepest 
descent [i^ from (I5T|) is it possible to prove that the asymptotic behavior of the purity for the QRW is oc l/y/Dt, 
this asymptotic regime is also shown (straight line) in figure 5. For the QLW (6 > 1) we can see that the asymptotic 
behavior looks like oc 1/t^, where the exponent ^ is a non-trivial function of Weierstrass' parameters h^A. In the 
inset of this figure we also show the long-tail exponent ^ versus the Weierstrass parameter A, here a transition in the 
behavior of ^ for h/A^ \ can be seen. From this inset it is also possible to see that if A ^ 1 the exponent ^ goes to 
0.5 corresponding to the QRW case. 



4- Quantum decoherence from a coherent preparation in \k) 



Now we will study the relaxation from a coherente preparation in the Fourier bases. In order to simplify this 
analysis we adopt here the following initial preparation for the system 



p(0) = -3- /" 'dfc /" " dk' \k) {k'\ , (32) 
kc Jo Jo 

where kc < n characterizes the initial uniform probability distribution in Fourier space. Therefore the mean- value of 
the Fourier modes appearing in the initial preparation is kpkk{0)dk — kc/2. On the other hand, from (j32p we get 



1 



Pk.kM ^ {ki\pi0)\k2) ^ - dk dk'S{ki-k)S{k' -k2) (33) 

fcc Jo Jo 

= -ie(fcc-fci)e(fci)e(fcc-fc2)e(fc2), 

where Q (z) is the step function. As indicated before the uniform initial probability distribution will be invariant in 
time pfcfc(O) — -^'d {kc — k)Q (fc), and fulfills normalization: 



/ dkpkk{0) = ^f Q{kc-k)e{k)dk^l. 

J — TT "-c J — 7r 



Therefore the probability {l\ p{t) \l) will have a time dependent spreading, but with a well defined mean-velocity 
characterized by the Fourier value kc- Using p3p as the initial preparation of the system, we get for the time- 
dependent probability to be at lattice site I 

/•TT /"TT 

{l\p{t)\l} ^ / dkj dk2{l\k,){ki\p{t)\k2}{k2\l) 



1 

2^ 



dki / dk2 ex.p [i {ki - k2) I] {ki\ p{t) \k2) 

— 7r J — 7r 



— — / dki / dk2exp[i{ki - k2)l]exp{I' {ki,k2,b,A)t) . (34) 
27rKc Jo Jo 
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FIG. 5: Log-log plot of the quantum purity Tr [p^] as a function of dimensionless time t when the system is prepared in 
the pure state p(0) = |Zo) {Zo|- The purity function is plotted for several values of Weierstrass' rate b/A. The purity for the 
QRW (straight line) shows the predicted long-time decoherent decay oc l/\/t, see ((31]). The long-time asymptotic behavior for 
the QLW shows an asymptotic behavior like oc where the exponent ^ depends on the Weierstrass' parameters fe, .4 in a 
non-trivial way, see inset. 



Note that at i = the initial distribution was delocalized in all the lattice according to: 

(^1 P(0) 10 = 7^ [ f ^ dk2 exp [z (fci - k2) I] (35) 

-iT^ric Jo JQ 

= ^ ^ (1 - coskj) , Ig{-oo,+oo). 

In figure 6 we have plotted the probability {l\p{t) \l) from as a function of the position I for different values of 
time. The non-diffusive behavior as well as the reentrance phenomenon of the (driven) profile probability {l\ p(t) \l) 
for the QLW can clearly be seen. What we call reentrance phenomenon is just the cooperative result of the clustering 
that occurs for b > A and the quantum coherence. Classically the clustering is certain when < ^ < 1; only in the 
case p > 1 the walker ultimately returns to fill in any gaps in the set of sites visited and the clustering eventually 
disappears [ssj . 

The agreement between the driven profile and the prediction associated to the pseudo-momentum is clearly visible 
from figure 6. Note that for the coherent preparation ()32|) . the thermal mean value of the pseudo-momentum will be 
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/ 

FIG. 6: Plot of the probability pi,i{t) = {l\p{t)\l) as a function of lattice site I for several values of dimensionless time t, 
Weierstrass' rate b/A, frequency rate r and for rc — ^, when the system is prepared in a Fourier coherent state like in (|32[) . 
The straight line corresponds to the profile at t = as is given in ([35} for kc = tt/4. figures (a) and (b) correspond to the QRW 
(fe = 1) for two values of rate r, the profile moves to the right with a thermal mean velocity as predicted in (|36|l . In figure (c) 
it is possible to see the " reentrance" of the probability to small values of I for increasing time {t — 20) , this phenomenon is due 
to the clustering of the visited sites and occurs above the critical value b/A > 1. 



different from zero 

(p(t)) = Tr[pp{t)]= I dk I dfci(fc|p|fci)(fci|p(i)|fc) 

/TT 
dki pfc^ (fci - k) {ki\p{t) \k) 
-TT 



1 



— / d/c pfe 

mn f A-l 



hkc \ A J ^ A 



J2 ±[1- cos {k,b-)], (36) 



where we have used ((26)) . Interestingly taking 6 = 1 in ((36)) we get that for the coherent initial preparation (|32t the 
QRW has a thermal mean velocity characterized by (p(O)qrw /"^ — W' ^'^ ~ (^c)] .Therefore we may conclude 
that for the initial preparation p2l) , the driven profile of the QLW moves faster than the corresponding profile for the 
QRW, as can be seen from the following inequality: 
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^ ^ n— 

Note from ([36|) that even when the eigenvalue pfc, see (f26| . is not defined for b > A, the double Fourier integral 
restores the regularity, i.e., the thermal mean-value of the pseudo-momentum (p(t)) is well defined for any value of 
the Weierstrass parameters b, A. Of course this is nothing more than the fact that by using the initial preparation 
([32]) the probabilistic profile move to the right with a finite velocity. In the case ^ ^ 1 the value of the thermal mean 
velocity goes to the one corresponding to the NN model. 



III. DISCUSSIONS 



In this work we have introduced the Weierstrass shift operators to build up the quantum master equation (also 
called Born-Markov equation) governing the quantum Levy walk process. We have started from the microscopic 
dynamics of a free particle, in a lattice, in interaction with a thermal phonon bath. The coupling with the heat bath 
was chosen proportional to shift operators (see (IC1[) and (IC2[) ). but any type other coupling can be worked out using 
the recipes that we have write down in appendices A, B, C. Then the Kossakowski-Lindblad generator of the quantum 
semigroup, which is the starting point to analyze the positivity condition of the structure matrix, was built up. For 
the present model it was not necessary to apply Davies's formalism (random phase approximation) to calculate the 
generator because Weierstrass' shift operators, in the Heiscnberg representation, are constant in time. This results 
in that the complete positivity condition of the generator was assured. The quantum Levy walk Hamiltonian shows 
for 6 > ^ that the DOS has a complex null-set inside the Brillouin zone, because the continuous eigenenergy £k{b,A) 
is non-differentiable when b/A > 1. We proved that several physical objects show the signature of this singular 
characteristic. A rigorous calculus based on the scaling properties of the eigenenergy (jS]) leads to the conclusion that 
the record of £k{b,A) as a function of k shows a critical behavior for values b > A. In fact, the eigenenergy £k{b,A) 
is self-affine with a (Box-counting) fractal dimension = 2 — /i, for < /i = \nA/\nb < 1. 

The analytical solution of the quantum master equation has been found in the Fourier basis, see (j22p . In particular 
the evolution equation for the second moment (HH), and its solution (pn|. has been calculated. Thus we prove that 
the thermal second moment is not divergent under the restriction b < A, which is different when compared with the 
classical Weierstrass universality b^ < A. The present quantum master equation can also be used to get information 
about higher moments and correlations of the quantum Levy walk process. The quantum decoherence has been 
characterized analyzing the probability to be in a lattice site considering different initial preparations for the system, 
i.e., we have worked out two different preparations for the reduced density matrix: a pure state in the Wannier basis, 
and a coherent preparation in the reciprocal Fourier basis. We have also examined the interference correlation function 
associated to the localization measurement process. In the present paper we work in a discrete infinite dimensional 
Hilbert space H with lattice parameter e = 1. The important issue about discrete against continuum dissipative 
model can also be tackled, for the present scale free Hamiltonian, by introducing the small lattice parameter limit 
e — > and using the Wigner transformation in the quantum master equation ^ (see end of appendix B for details). 

The continuous eigenvalues of the scale free Hamiltonian ([1]), and the effective Hamiltonian dTU]) have been studied 
as a function of the Fourier number k (the tight-binding case corresponds to take 6 = 1 in ([5])). The quantum 
Levy walk Hamiltonian (i?s) for different values of Weierstrass' parameters b and A has been plotted showing the 
non-differentiable structure of its spectrum £k{b,A) when the rate b/A is larger than 1, this was also quoted in 
connection to the critical behavior of the DOS ([T]). We show that the continuous eigenvalue of the effective Hamiltonian 
{Hs — hojc a'^a) has fluctuations of larger amplitude for increasing energy rate rg = huc/^, see figure 2. 

For the transport of a quantum Levy particle the non-trivial effective Hamiltonian can be characterized in terms of 
the frequency rate Tc = uJc/D, this dimensionless number compares the Caldeira-Legget upper frequency fl^ against 
the diffusion coefficient D = T:a/2ph, being a the dissipative parameter and /3 the inverse of the temperature of the 
bath (see lC7p . When the system is prepared in a pure Wannier state, the probability to be at distance from the 
initial position Iq, i.e., P{Al,t), was studied as a function of Al for several values of time t and for several values of 
the frequency rates r = ^/D, r^,. The larger the frequency rate r is the longer the quantum coherence persists before 
the profile starts to spread. The opposite case r — >• corresponds to the high temperature limit. We have shown that 
the profile P(Al,t) spreads in a non-diffusive way when & ^ 1, i.e., characterizing a Levy-like behavior. The larger 
the rate Tc is the more the time oscillations of P{Al = 0,t) there are. We have also shown that the profile shows an 
important trimodality, this behavior is intensified by increasing Weierstrass' rate b/A above the critical value b/A — 1. 
The localized correlation x{t) has been studied as a function of time t when the Levy particle is prepared in the pure 
state p{0) = \lo) {lo\. The function x(i) has been plotted for different values of Weierstrass' parameters b,A and 



18 



frequency rate Tc, we have shown that the asymptotic decay of x{t) does not have the ID classical oc 1/v^ diffusive 
behavior, on the contrary the localized correlation function x{t) shows a long-time coherent persistence, see figure 4. 
The quantum purity Tr [p^] has been studied as a function of time t for several values of Weierstrass' parameters 6, A. 
The purity for the quantum random walk (6 — 1) shows the predicted long-time decay oc for the NN case (|3ip . 

Nevertheless, for the quantum Levy particle the purity has a long-time asymptotic behavior that looks like oc l/<^, 
where the exponent ^ depends on the Weierstrass parameters b/A in a non-trivial way, see inset of figure 5. 

When the system is prepared in a coherent state as in (I32[) . the probability (/| p(t) \l) for the quantum Levy particle 
has been studied as a function of I for several values of time t, Weierstrass' parameters 6, A and frequency rate r. 
We have checked numerically that for the coherent preparation (l32)l . with kc = 7r/4, the profile moves to the right 
with a thermal mean- velocity as predicted in (1361) . We have also shown that the profile has an interesting reentrance 
behavior. In figure 6c the reentrance of the probability to small values of I for increasing time is clearly visible, this 
phenomenon is the result of the cooperative phenomena of clustering (for the critical value b > A) and quantum 
coherence. The trimodality that occurs from a pure initial preparation like p(0) = |Zo) (^oL also is shown to occur 
when there is clustering in the walks. 

In the present framework it is also possible to analyze the quantum jump picture (jASp which is of great value to 
measure the dissipative decoherence [H, . This last analysis can be done by studying the fluctuation superoperator 
F [•] in the Fourier basis, i.e., (|C7p . We have prove in the appendix C, that the superoperator F [•] can be handled 
in terms of the elements of the operators a,a^, which for b ^ 1 have a non-trivial Fourier structure resembling the 
characteristic function of the classical Weierstrass random walk. Further works along these lines are in progress. 

In conclusion, the open quantum Levy lattice model introduced here for transport in a infinite dimensional Hilbert 
space, is of remarkable usefulness for calculating analytically many interesting measures of irreversibility. The impor- 
tant issue of decoherence in a dissipative long range model was tackled analytically. The analysis of other irreversible 
measures like entropy in the context of the coherence- vector formulation, etc. |2j|, can also be worked out in the 
present framework and is in progress. 
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Quantum dynamic semigroups are the generalization of Markov semigroups for non-commutative algebras [ig . |2J] . 
In the Markovian approximation, Kossakowski and Lindblad established the form of the Quantum Master Equation 
(QME) in order that the evolution of the system of interest shall correspond to a quantum dynamic semigroup 
(these semigroups are also called "completely positive semigroups" CPS). In the structural theorem, Lindblad [26[ 
established that the generator of a CPS, acting on the reduced density matrix of the system, p, has the form: 
L [•] = — ■|-[_ffe//, •] -I- ^X^Q [^a'j^a] + [^ai'V^cJ], wheic Va are bounded operators. If the system has a finite 
number of degrees of freedom, the expression for the generator L [•] can alternatively be written in a different way. 

Consider for example the algebra of the N x N complex matrices, then assume that the set of operators {Va}^^i^ 
fulfilling Tr[Vj^VQ'] = Sa,a', is a basis in that space (i.e.: the basis is orthonormal with respect to the scalar product: 
< A\B >= Tr[A^B]). In term of this basis, the generator of the CPS is [25,] 



where the matrix of elements aa-y is positive-definite. This generator is written in the Schrodinger representation and 
acts on the density matrix of the system of interest. In the Heisenberg representation the dual generator L* [•] , defined 
as Tr[L*[A]p] = Tr[AL[p]\ acts on any physical observable A (Hermitian operator). Now we define the superoperator 



Appendix A: Quantum dynamic semigroups revisited 




(Al) 




(A2) 



and we consider its dual F*[*] evaluated in the identity operator 1, i.e., F*[l] — ^ 
(|A1|) can be written in the compact form 



a, 7—1 



aa-yV^jVa ■ Then the generator 



L[.] = -^[H,ff,>]+F[.]~^{F*[l],.}^, 



(A3) 
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The operator if can be regarded as a dissipative operator, and the fluctuating superoperator. Equations 



[2]) and (IA3|) allow us to define the structure matrix [oq,^]. This matrix contain the information concerning the 
relaxation times of the dynamic system in contact with a thermal bath. We say that a generator with the structure 
()A3|) . with an Hermitian matrix [oa^] has the form of a Kossakowski-Lindblad generator (KL) |28]. We note that the 
corresponding semigroup is completely positive if and only if [aa-y] is a positive-definite matrix, and it is equivalent 
to say that the generator L [•] fulfills the structural theorem. Alternatively, we will say that when [aaj] > the 
generator is a well defined KL generator 44, 45]. 
The formal solution of the QME, p ~ L[p], is 

Pit) = exp + F[.] - i {F*[l],.}^)t^ piO). (A4) 

This expression can be put in the form [l^ 



x{S{t - t,n)F[*]S{tra - tm-l) ■ ■ ■ F[*]S {ti)} p{0) . 



In this way the dynamics of the system can be interpreted as if it were composed by quantum jumps (associated to 
the superoperator F[u]) and in between them there is a smooth non- unitary evolution determined by 

Sit)p = exp .] - I •}+) t| P ^ N(i)pNt(i), 

where N(i) = exp(— [j^Heff + t) characterizes decay. This representation is very suitable for describing the 

decoherence of the off-diagonal elements of the density matrix [2^ |4^ . 



1. The quantum master equation and the second order approximation 

It is known that the QME arising from second order perturbation theory has, in general, the KL form (|A3p . To 
see this assume that the total Hamiltonian is of the form: Ht — Hs + Hb + OHj, and that the system S interacts 
with a equilibrium thermal bath B through the term OHj (9 is the coupling intensity) , also we assume that the initial 
condition for the total density matrix can be written in the form pt(0) = p{0) (g) pg, where pg is the equilibrium 
density matrix of the bath. Now consider the Liouville equation for the total density matrix and trace out the bath 
variables, keeping only up to the second order 0{9^). This procedure gives a QME for the reduced density matrix 
of the system p having a KL form where the generator is defined through an effective Hamiltonian Heff and the 
superoperator F [•] [23] , 

H,j^=Hs~i-j^ dTTrB{[Hi,Hi{-T)]p%), (A6) 



F W)] - f - j dr TrB (Hj p{t) ® p% Hi{-t) + Hi{~t) pit) ® p^ Hi) . (A7) 

Where Hj{—t) = e-«'^(^fs+-H^B)/'i//^ gi-r(_ffs+_f/B)/;i_ remark that this structure for the generator is independent of 
any particular system S under consideration; it is also valid for finite or infinite dimensional Hilbert spaces. Now we 
consider the interaction Hamiltonian to be characterized by the direct product of operators: 

Hi = Y,yp®B0, (A8) 



Then using explicitly the Hermitian condition of Hj, H^ff and the form of the superoperator F [•], we can write 



n2 roo 

H^ff ^Hs-i-Y, dr [xo.p{-T)V^Vp{-T) - x:^(-r)y;(-r)K. 

aP 



(A9) 
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(AlO) 



Here we have introduced the correlation functions of the thermal bath: 

Xaphr) = Ti'B {p'b Bi Bpi-r)) , (All) 

where 

Ba{—T) = cxp{—iTHB/h)BaC:iq){iTHB/fi), 
V^(— r) = exp{—iTHs /h)Vi3 enpiirHs /h). 

Because the thermal bath is stationary the correlation function fulfills the symmetry condition: Xq/3 ("''') — Xpai''')- 
The KL form (IA3|) allow us to analyze its possible positivity. We do this because, even when the QME up to the 
second order approximation can be written in a KL form, it is not possible to assure that the semigroup will be 
completely positive. 

In the case of working with a finite dimensional Hilbert space the analysis of the structure [oa-y] allows us to 
introduce a necessary condition on the Hamiltonian Hj in order to arrive to a well defined KL. Assuming, that the 

n 

interaction Hamiltonian can be written (in any particular basis) in the form: Hj = Bp with n < N'^ — 1. 

The set {V^}^^]^ must be closed in the Heisenberg representation, i.e.: 

rn 

Va(-T) = C0^{-T)Vy with m<n; (A12) 

7=1 

otherwise the matrix [ua^] will not be positive-definite |28(1. If the KL generator were not a genuine CPS we ought to 
apply some random phase approximation (Davies' devices [27j). 



Appendix B: On the Weierstrass shift operators 

Consider the product of two Weierstrass' shift operators, using we write (for A > l,fe(integer) > 1) 

_^ \ 2 CXD oo oo oo 



ni—0 1— — 00 n2—0 I4 — — 00 



Then, in Wannier's representation we get the off-diagonal elements 

^ ^ ni,n2— ^,/3,/4 — — 00 



A-V ' °° 



A J ^ ^ A 

^ ni,n2=0 £3 = — 00 

E ^».+»2 ' iA>i,b>i). 

ni,n2— 

In a similar way it is simple to show that: 

{h\a^a\l,)=f^y E '-^"^^^^ (-4>M>1), (Bl) 

therefore {li\ aa^ IZ2) = {li\ a'I'a IZ2), telling that [a'' ,a] =0. ■ 

Note that if 6 = 1 we explicitly get from (|B1I) that a^a — 1. Assuming that b — integer, we can write for the 
diagonal elements 



^<L (.4>1,6>1), 



A I ^-^ Ani+n2 , 

^ «i,n2=0-^ I 1, (yl > 1,6 = 1) 
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Alternatively, in the Fourier representation we can write: 



{ki\a^a 1^2) = ^ (fci - ^2) 



A-1 
A 

A-l 
A 



exp[-?A;i (6"2 - 



E 

Tl2,^ll— 

' ^ COSfci _ 



n2 



(B2) 



So if 6 7^ 1 (i.e., for the QLW) a^a is not the identity operator, only in the NN case {b — 1) we get a^a — 1. From all 
these results and considering the structure of the Hamiltonian Hs we may conclude also that [Hs , a^] = [Hs ,a] — 0, 
telling us that Weierstrass' shift operators, in the Heisenberg representation, are constant in time. 

Here we calculate the matrix elements of [q, Hs] for the general case b ^ I. In the Wannier basis we obtain 

{h\[<l,Hs]\l2) = ^ai|q|0a|i/5|/2>"(;i|i?s|0('|q|'2) 
I 

= J2lSi,j{l\Hs\h)~hSLiAh\Hs\l) 
I 

2A [^A-^ 



n f A-i 



/2.-b" 



2 \ A J \A 
Note that when 6 = 1 we get [q_,Hs] — 5^/2 (a — a^), because 

n ( A-\ 



(B3) 



A 



^E (;^) ((^ii-is -1 - <5ii-i2,i)j 



§ai|(«-«0l^2) 



Therefore the operator p — [a — a^) = (e"*^^' — e'^^') can be associated to a "discrete" momentum operator 
in a lattice with scaling parameter e. Also it is trivial to see that [p,_ff5] = as expected for a free particle model. 
As a matter of fact, taking the limit of the lattice parameter going to zero, i.e., 

lim = limcosh(e9;) ^ (1 H — e9;), 

e— ^0 2 e— >0 2 



lim 

e->0 



(a -at) 



— lim sinh(ec};) — > —{edi), 



we recover the usual commutation relation [q, p] = ^ ^ 



(a+a+) 



ihi [13]. 



Appendix C: The QME for the density matrix of the QLW 



Here we find the QME for our quantum Levy model. To write the QME we have to calculate the superoperator 
and the effective Hamiltonian H^ff, both objects are given in (|A7I) and (jA6l) respectively. We assume that the 
interaction Hamiltonian (jASp is written in term of two system operators 



Vi = hVa = v^, r > 0, 



(Cl) 



and two baths operators (infinity set of thermal harmonic oscillators [47[) 



Bl 



(C2) 
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which fulfill Xq^(— t) = Ttb {p% B^{—t)^ = Sa^Xaai—T). So the correlation functions of the bath are character- 
ized by: 



Xi(-t) = ^\vkfexp{-iuJkT){n{uj) + l), 
k 

X2{-t) = ^\vkf exp{iuJkT)n{uj), 



(C3) 
(C4) 



where n{uj) — [exp{haj/kBT) — 1] ^. Using that Vq(— r) = V'ct(0),a = 1,2, i.e., they are constant in time, we can 
write from (jAlOp the fluctuating superoperator in the form 



(hrf / dr [(xi(-t) + xli-r)) a.a^ + (x2(-r) + xli-r)) • a] 



and from (IA9I) the effective Hamiltonian as 



(C5) 



n2 roo 

Heff-Hs = -'^T.J^ dTix^„i-T)~x:c.i-r))V^V^ 



(C6) 



= _i_(ftr)' 



dr (xi(-t) - Xi(-t)) a)a + (x2(-t) - X2(-^)) 







From (jCSp and (jC4p the Fourier transform of the thermal bath-correlations are 



hi{uj) 

/l2(w) 



dr exp(— iw)xi(— t) = 27rg(— w) oj) -f 1] 



— oo 

+ 00 



dr exp(— jcj)x2(— t) = 27r5(cj)ri(a;), 



where g{Lo) — X]fc ("^ ~ '^fc) the spectral function of the phonon baths [13l |. In general the half-Fourier 

transform, that appear in (|C5p and lC6l) . can be written in terms of ha{uj), i.e.,: 

(irexp(--ia;)xaQ(-r) = — hi — ^ — , 

where Sa{i^) — I^oo'^^ ha{u){u — oj)"^ is its Hilbert transform. Note that h^iuj) € TZ,,, and so Sq(w) G 7?.e, 
because the thermal bath correlation function is stationary Xaai~T) — Xaai^)- Then we only need to evaluate the 
integrals 

poo 

dT(Xaa(-T)+xL(-T)) = h^{0), 
dr (Xaa(-r) - Xaa(-T-)) = iSa(O). 

The c-numbers {ha{0), Sa(0)} are given in terms of the Fourier representation of the bath correlation function at zero 
frequency. Modeling the coupling constant in the Ohmic approximation g{uj) = goj if < w < Wc, 5 > (where 
ujc is larger than any characteristic system's frequency [l^) and taking into account that wn(cLi)|^^p = {j3h)~'^ , we 
get hi{uj ^ 0) = h2{'-o ^ 0) = 2'Kg{j3h)~^ where /? = l/ksT is the inverse of the temperature of the bath, and 
si(lo — > 0) = S2(w — > 0) = —2gLbc- Then 



= 0^T^2TTg {pny^ [a • a'' + a'' • a] 

Tra r t t 1 n 
= ^ [a» a' + a' • aj , a > 0, 



(C7) 
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where a = 4:9 T g is the dissipative constant. For the effective Hamiltonian we get 

Heff = Hs - — ihVf {siiO)a^a + S2{0)aa^) 

= Hs- 2g(JjcO'^r^ha^a ^ Hs - hwcO^a, > 0, (C8) 

where tOc = 2gujc0^T^ is an upper bound frequency. With these expressions for and H^ff we can write, using 
(|X9)) . (|MOt and the QME in the form 

p = :^[H,.jf,p] + F[p]-^{F*[llp}^ (C9) 

" T t-^'^//' ^] + ^ (Sapa'f - a'^ap - pa'^a) + ^ {2a'' pa ~ aa'' p - paa^) . 

Now we exphcitly calculate the evolution equation for the elements of the density matrix p. Using the Fourier basis 
in the we get 

d ~~ % IT (y. 

^(fci|p|fc2) - —{k^\[H,ff,p\\k2) + — {k^\{apa^ + a^pa)\k2) (CIO) 

\ {aa'' p + paa'') jfca) , 

where H^ff = Hs — hucCi^a, and we have used that [a, a^ — 0. Note that 

(fcilaaVlfe) = J dk{ki\aa'<\k){k\p\k2) (CU) 

^ ^ n2,ni— 

^-iV ^ cosA:i(6"^ , ,, , 

2. j;^, (fell P 1^2), 



_4 / ^-^ J[ni+n2 



and 



(fci|paat|fc2) = r dk{k,\p\k)Sik2-k)f^y £ (C12) 
'^-ly cosfc2(6"--b"0 , 

n2,ni—0 



where we have used (jB2p . Using that 



{h\a\k2)^(^)±'^^Sik.-k2). 



^ A J A" 

we can also write 

{ki\apa^\k2) ^ I dk I dk3{ki\a\k){k\p\k3){k3\a^\k2) (C13) 



— 7T <J — TT 



^— 1\ ^ expzfcife^^ exp — 7^26^ I I, . 
^) 2^—4^^ {ki\p\k2) 



and 
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Therefore putting (IC11IC12IC13IC14P in (p^ we get the final result (for 6 > 1,^ > 1) 

Jt{ki\p\k2) = ^(E,,-E,J(fci|p|fc2) 



na ( A—1 



2 oo 



E 7^[2cos(fci6"-fc2n 



2hl3\ A J "^^A 

^ ' n.rn—{) 

■ cos kx (&" - 6") - cos k2 (&" - 6")] (fci \p\k2), 



where E^^ is the effective eigenenergy: 



/ 4-1 1 \ 
E., ^ l^(^l--^^^_cos(5"/ci)j 

^ cosfci 

^ __ . 

n2,ni— 



(CIS) 



(C16) 



Appendix D: The second moment of the QLW 



In the present appendix we explicitly calculate the quantum thermal mean- value (c^{t)^, thus 

/>7r />7r 

(q^(i)) = Tr [qV(i)] = E / / ^^^1*2 (/| |fci) {k^\ p{t) jA^) (fca \l) , 



(Dl) 



where (fc2| /3(i) l^i) is known form our general solution Eg. ([2^ . Using the pure state as the initial condition for 
the density matrix we get 



(q2(0) = f dki f dk2{l\q'\li){li\k,){ki\p{t)\k2){k2\l) 

- ^ E f f ?'<5z.iie^'^'i (All p(0 |fe2) 6-*'=^' 

^ /'TT /'TT 

rffci(ifc2 ;2g«(fei-fe2)(i-io) exp(J'(fci,/c2,6,^)i)- 



^ /'TT /•TT 



(27r)2 E 



Introducing the change of variable I — Iq ^ using the "functional" properties of the Dirac-delta [48 1 

^ oo 

ZTT ^ — ^ 

/ — — oo 

/(fc)<5(")(fc-fco)dfc = (-ir/(")(fco), 

and after integrating by part we finally get the result 



(D2) 



<q'(i)) = ^ E / / (; + Zo)'e'('=i exp(^(fci,fc2,&,^)i) 



— TT J —TT 



(fT{k,k2,b,A) 



«7r 



dk^ 

( dJ^ {k,k2,b,A) 
\ dk 



dT{k,k2,b,A) 



dk 



k=k2 



k=k2 



72 

27r 
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where we have used that exp {T {k2, k2,b, A) t) = 1. By introducing the exphcit expression (|23)) . for {k,k2,b,A) , 
and noting that we have assumed b = integer, we finally obtain (|20p . 

As we have already pointed out when working on the differential equation for the second moment, see (|f4|) and 
([T^. here we have explicitly shown that only if 6 < ^ we get (q^(t)) ^ oo, this conclusion is quite different when 
compared to the classical counterpart of a Levy flight, where a divergent second moment appears if 6^ > A. The 
reason is that for the QLW the second moment (q^(t)) cannot just be calculated from the second derivative of the 
Fourier transform of the space-probability distribution [H, 0] ■ In quantum mechanics the thermal second moment is 
given through a trace operation, and {ki \ p(t) \ k2) takes into account all the coherence phenomena. 



Appendix E: Density of relaxation for a classical homogeneous walk and its localized probability 



Consider a classical one dimensional homogeneous random walk, its master equation will be 

(s, t\ 0,0) =D [P (s + 1, t\ 0,0) + P (s - 1, t\ 0, 0) - 2P (s, t\ 0, 0)1 , P (s, 01 0, 0) = Ss o- 

at 

Introducing the Fourier transform of the probability A (fc, t) = X]^-oo e^^^P (s, t\ 0, 0), we get the solution A (fc, t) — 
exp{2Dt [cosfc — 1]) ■ Therefore the relaxation form the initial condition (classical localized probability function) will 
be 

Po{t) = — \{k,t)dk= e-^'p{j)dr, 7>0. (El) 

Where ^(7) is the (probability of) relaxation density: 

p(7) = — / 5 {-f- 2D [1- cos k]) dk ^ -UjD-j^)~^'^ , 
2n n 

in total analogy to the density of states of the ID tight-binding model [3, '9']. From the expression ^(7) and Eq. (IE1[) 
it is simple, by using the Laplace transform, to show that asymptotically for t — )■ 00 

Po(i)^^, (E2) 

this is the typical behavior for the ID localized probability in a classical diffusive regime. 



